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Rapidly quenched Kosterlitz-Thouless (KT) superfluid transitions are studied by solving the Fokker-Planck 
equation for the vortex-pair dynamics in conjunction with the KT recursion relations. Power-law decays 
of the vortex density at long times are found, and the results are in agreement with a scaling proposal 
made by Minnhagen and co-workers for the dynamical critical exponent. The superfluid density is strongly 
depressed after a quench, with the subsequent recovery being logarithmically slow for starting temperatures 
near Ta't- No evidence is found of vortices being "created" in a rapid quench, there is only decay of the 
existing thermal vortex pairs. 
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There is considerable current interest in the properties of 
quenched superfluid transitions, since they may have rele- 
vance to cosmic-string phase transitions in the rapidly-cooling 
early universe [p. Experiments [||] have been carried out in 
superfluid 4 He attempting to observe vortex structures "cre- 
ated" in rapid pressure quenches through the transition, but no 
vortices have been observed, which is a puzzle [||. 

To try to understand these quenched transitions we have 
constructed [Q] an analytic formulation of the quenched 
Kosterlitz-Thouless (KT) superfluid transition in two dimen- 
sions (2D). Although not the three-dimensional (3D) realm 
of the experiments, the simplicity of the KT vortex renor- 
malization scheme allows insights into the vortex dynamics 
that should be useful in extending the calculations to 3D 
Previous studies of quenched 2D transitions have included 
computer simulations [|6j ^, |8|] and dynamic scaling theories 
||^, |T(], [HI], but there is still considerable controversy over the 
physics involved. The various dynamic scaling theories give 
quite different predictions for the temperature dependence of 
the dynamic scaling exponent z, and Bray and co-workers 
0] have even suggested that there may be violations of dy- 
namic scaling. Unfortunately the computer simulations are 
not able to resolve the questions in the field, since it has re- 
cently been shown H that different boundary conditions in 
the simulations give different results for the dynamics. 

We consider a thin superfluid helium film on a flat substrate, 
which is coupled to a thermal reservoir at temperature T. The 
vortex pairs in the film are characterized by the distribution 
function T(r,t), which is the density of pairs of separation 
between r and r+dr. It is determined from the 2D Fokker- 
Planck equation [|l2|], 
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where r 1 = r/a Q with a the vortex core radius, U is the 
pair interaction energy, and D is the diffusion coefficient char- 
acterizing the mutual friction drag force on the vortex cores 
of a pair. It can be seen that the time in this equation is 



scaled by r D = a 2 /2D, the diffusion time of the smallest 
pairs of separation a a . In common with the scaling theo- 
ries and simulations, we assume that any renormalization of 
D can be neglected. At equilibrium where ^7 = the 
solution of Eq. (fil) is just the usual Boltzmann distribution, 
T = a~ 4 exp (-{U + 2E c )/k B T) where 2E C is the core en- 
ergy of a pair. 

The vortex interaction energy U in Eq. ([]]) is determined 
using the Kosterlitz-Thouless (KT) vortex renormalization 
methods [O]. The Kosterlitz recursion relations can be writ- 
ten in the form 
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where K = h 2 a s /m 2 k B T is the dimensionless areal super- 
fluid density. We take the core energy to be the Villain-model 
value 2E c /k B T — tt 2 K , where K a is the starting value of 
K at r = a Q . In thermal equilibrium Eqs. (]l]-^|) lead to the 
well-known result that the superfluid density has a universal 
jump to zero just above the transition temperature Tkt [0- 
To study the quenched transition [[ij the film is first equi- 
librated with the heat bath at temperature T, generating an 
equilibrium distribution of vortex pairs. The top curve in 
Fig. |] at t = is the distribution at Tkt', it varies asymp- 
totically as (r I a Q )~ 2 ' lxK . The temperature of the heat bath is 
then reduced abruptly to a low temperature, 0. 1 Tkt, where in 
equilibrium the vortex density is over 20 orders of magnitude 
smaller. However, due to the Kelvin circulation theorem the 
pairs cannot suddenly disappear, since the only way they can 
be extracted by the heat bath is if the plus-minus pairs annihi- 
late at r = a , where the remaining core energy is converted 
to phonons. This can take quite a long time to occur, however, 
since the pairs only slowly lose kinetic energy to the diffusive 
frictional force of the heat bath as they move together towards 
annihilation [ 15 1. 
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To solve for the time dependence of the distribution func- 
tion after the quench we combine Eqs. (Q) and (||), which be- 
come 



they screen more effectively due to the lack of thermal fluctu- 
ations. 
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We solve this in conjunction with Eq. (fy using the method 
of lines with a third-order Runge-Kutta technique on a fi- 
nite domain, to a maximum pair separation R/a = e 10 = 
2.2xl0 4 . The boundary conditions used are that the flux of 
pairs J = —2D(r^ + 27rKT) be continuous across the 
boundary r = R, and that T drop abruptly to zero at r = a Q , 
which is equivalent to putting a delta-function sink term at a 
on the right-hand side of Eq. (Q) to generate the annihilation 
at that point. 
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FIG. 1: Time dependence (in units of the diffusion time r Q ) of 
the pair distribution function T (in units a^ 4 ) for an instanta- 
neous quench from Tkt to 0.1 Tkt- 

Figure |l] shows the time dependence of the distribution 
function for an instantaneous quench to 0.1 Tkt from a start- 
ing temperature barely below Tkt (from K a = 0.747853, 
compared to the critical value 0.747852 at Tkt)- As expected 
for the annihilation process the smallest pairs decay away first, 
in the classic pattern of a phase-ordering transition where at 
long times the largest scales become more dominant. By in- 
tegrating the distribution function over d 2 r the pair density is 
obtained, shown as the top curve in Fig. ^a. At long times the 
density decreases nearly as t , in agreement with the scaling 
theories [l(| [□} at Tkt- We find, however, that the slope 
only approaches -1 logarithmically slowly, being -1.12 at t = 
100 and -1.09 at t = 2x 10 4 , and more detailed fits show that 
the density is varying as (tint) -1 . This is a consequence of 
the slow recovery of the superfluid density after the quench 
(bottom curve of Fig. ^b ), which has a logarithmic approach 
to its equilibrium value. The substantial drop in the superfluid 
fraction just after the quench (from an initial value 0.89) is due 
to the increased polarizability of the pairs at low temperature. 
Their density is unchanged immediately after the quench, but 
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FIG. 2: a) Decay of the pair density (in units a~ 2 ) and b) re- 
covery of the superfluid fraction (evaluated at r = R) in time 
after a quench to 0. 1 Tkt from different starting temperatures 

Quenches from starting temperatures below Tkt also dis- 
play power-law variation of the vortex density, since the KT 
transition is a line of critical points. The plots in Fig. [l]show 
this for a series of initial temperatures between 0.8 Tkt and 
Tkt- At the lower starting temperatures the superfluid den- 
sity recovers more quickly, and the vortex decay then varies 
accurately with the form t~ z / 2 , where the constant z agrees 
with the scaling prediction of Minnhagen and co-workers for 
the dynamic exponent [O], 
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Figure |3| shows the values of z extracted from the slopes of 
the curves in Fig. which are compared with z sca ie, with a 
dynamic exponent predicted in Ref. [J12], and with simulation 
results. It has recently been shown [p] that the simulations 
of the vortex dynamics can be misleading, since different re- 
sults for z are obtained depending on the boundary conditions 
used: periodic boundary conditions (PBC) yield ^ only the 
spin-wave result z = 2 at all temperatures, while fluctuating- 
twist boundary conditions (FTBC) give z sca [ e 081. Our results 
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in agreement with z sca ie (except near Tkt where the logarith- 
mic corrections enter) provide further support that this is the 
correct exponent characterizing the vortex dynamics. Unlike 
the simulations our formulation does not contain spin waves, 
and this allows a direct probe of the vortices, without their 
behavior being obscured by the slower dynamics of the spin 
waves. 
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FIG. 3: Dynamic exponent z calculated from the slope of the 
curves in Fig. ||a, compared with various scaling theories. The 
PBC simulations are from Ref. 0, the FTBC from the relax- 
ational dynamics of Ref. [||]. 

Bray and co-workers JlO[ | have predicted the appearance of 
a logarithmic term in the decaying vortex density for quenches 
starting above Tkt, but with a form hxt/t that differs from 
ours. In their calculation they assume, however, that the su- 
perfluid density is a constant in time. As seen in Fig. ||b this 
may be a reasonable approximation for quenches starting well 
below Tkt, but it appears to fail near Tkt because of the 
logarithmic recovery of the superfluid density. The log term 
in our vortex density arises entirely from this variation. In 
Ref. [0 the appearance of the logarithm is ascribed solely to 
the presence of free vortices, and this is cited as a violation of 
dynamic scaling. This is not the case in our calculation, where 
no free vortices at all are present. It appears that the develop- 
ment of the log term only depends on a sufficient density of 
pairs of large separation; for quenches starting from T/Tkt = 
0.98 we find that the deviations from the power-law decay are 
already strong, but not quite the logarithmic behavior found 
right at Tkt- 

Figure ^ illustrates the approach to equilibrium following 
an instantaneous quench from Tkt to 0.75 Tkt- In the dis- 
tribution function the smallest pairs quickly come to thermal 
equilibrium (the dashed line) at the new temperature, while 
the larger pairs only slowly readjust. The vortex density lev- 
els off near the equilibrium value, since it is dominated by the 
smaller pairs. The superfluid density, determined by the larger 
pairs, recovers more slowly, although because of the higher 
temperature its initial reduction is considerably smaller than 
that in Fig. §b. 



FIG. 4: a) Pair distribution function and b) vortex density and 
superfluid fraction (dotted curve) versus time for a quench 
from T KT to 0.75 T KT - 

From Fig. [|a it is clear that no vortices are being "cre- 
ated" in the instantaneous quenches from Tkt and below; 
there is only decay of the existing thermal vorticity. We have 
looked at the effect of the quench time on this decay, shown in 
Fig. U This shows the vortex density and superfluid fraction 
for the case of a linear ramp of the temperature from Tkt 
to 0.1 Tkt in a time tq. We assume for convenience that 
D is independent of temperature, although in a real helium 
film it would vary during such a finite-time quench JT^]. It 
can be seen in Fig. [|a that for a finite quench time the vor- 
tex density decreases more slowly than for the instantaneous 
quench (dashed curve), since the system spends more time at 
higher temperatures where the thermal density is higher. As 
the temperature falls and the smaller pairs that had been able 
to stay in quasi-equilibrium start to annihilate, the more slug- 
gish larger pairs begin to dominate, and the curves revert to 
the instantaneous-quench result at long times. The slight wig- 
gles in the curves are due to changes in the superfluid fraction, 
shown in Fig. ^p. The fraction begins to increase slightly as 
the density of the smaller quasi-equilibrium pairs decreases 
with the decreasing temperature, but finally the increasing po- 
larizability of the larger pairs overcomes this, driving down 
the superfluid fraction. Since there are now more pairs re- 
maining than in the instantaneous quench the fraction drops 
even below the value for that case, and then rejoins it at long 
times as the smaller pairs decay completely and the larger 
pairs dominate. 

The behavior of the vortex density in Fig. |^ is the opposite 
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FIG. 5: Vortex density and superfluid fraction versus time for 
quenches from Tkt to 0. 1 Tkt in time tq (in units of the 
diffusion time). 



of the Kibble-Zurek predictions of Ref. [Eh, which postulate 
that the faster the quench is performed the more vortices will 
be created. Since these theories do not take into account the 
thermal vortices it is unclear how any extra vortices would 
be formed. Of course, our calculations cannot be extended 
to starting temperatures above Tkt, since the KT recursion 
relations become invalid. We doubt, however, that the addi- 
tional thermal vorticity generated at higher temperatures will 
do anything but add smoothly to the decaying curves of Fig. 
|[ The 2D quench simulations from above Tkt ® show no 
evidence of any additional vortex creation. Even simulations 
of quenched 3D transitions show (in Fig. [j]of Ref. Jl7[]) only a 
smooth and monotonic decay of the thermal vorticity starting 
from above T c , with no extra vorticity being generated. 

In summary, an analytic formulation of 2D quenched su- 
perfluid transitions shows a power-law decay of the vortex 
density with time, in agreement with the predictions for the 
dynamic scaling exponent of Minnhagen and co-workers. 
Quenches starting near Tkt show an additional logarithmic 
variation resulting from the slow recovery of the superfluid 
density to its equilibrium value. The strong depression of 
the superfluid density following a quench has not been pre- 
viously been taken into account, though it plays a major role 
in the quench dynamics, as seen in Fig. |[ It will be important 
that any future simulations of quenched transitions (in both 
2D and 3D) include computations of the superfluid fraction 
(the helicity modulus) as well as the vortex density, and also 



study the effect of changing the boundary conditions. Experi- 
mental tests of these vortex dynamics should also be possible: 
measurements on Josephson-junction arrays have been able 
to determine z = 1.98± 0.03 at Tkt [ 1 8 1, and experiments on 
liquid-crystal films are able to visualize the pair-defect anni- 
hilation process following a quench [JfJ]]. 

This work is supported by the National Science Foundation, 
DMR 97-31523. 



REFERENCES 

[1] A. Gill and T. Kibble, J. Phys. A 29, 4289 (1996); W. 
Zurek, Phys. Rep. 276, 177 (1996). Extension of this the- 
ory to 2D superfluids is found in A. Yates and W. Zurek, 
Phys. Rev. Lett. 80, 5477 (1998). 

[2] M. Dodd et al, Phys. Rev. Lett. 81, 3703 (1998). 

[3] Very different experiments in 3 He have claimed results 
consistent with the Kibble-Zurek predictions of Ref. 1 : 
V. Ruutu et al, Phys. Rev. Lett. 80, 1465 (1998); C. 
Bauerle et al, Nature 382, 332 (1996). 

[4] H. Chu and G. A. Williams, Physica B 284-288, 115 
(2000). 

[5] G. A. Williams, J. Low Temp. Phys. 93, 1079 (1993); G. 

A. Williams, Phys. Rev. Lett. 82, 1201 (1999). 
[6] R. Loft and T. DeGrand, Phys. Rev. B 35, 8528 (1987); 

M. Mondello and N. Goldenfeld, Phys. Rev. A 42, 5865 

(1990); B. Yurke, A. Pargellis, T. Kovacs, and D. Huse, 

Phys. Rev. E 47, 1525 (1993). 
[7] H. Luo and B. Zheng, Mod. Phys. Lett. 11, 615 (1997). 
[8] B. Kim, P. Minnhagen, and P. Ollson, Phys. Rev. B 59, 

11506 (1999); L. Jensen, B. Kim, and P. Minnhagen, 

Phys. Rev. B 61, 15412 (2000). 
[9] A. Rutenberg and A. Bray, Phys. Rev. E 51, R1641 

(1995). 

[10] A. Bray, A. Briant, and D. Jervis, Phys. Rev. Lett. 84, 
1503 (2000). 

[11] P. Minnhagen, O. Westman, A. Jonsson, and P. Olsson, 

Phys. Rev. Lett. 74, 3672 (1995). 
[12] V. Ambegaokar, B. Halperin, D. Nelson, and E. Siggia, 

Phys. Rev. B 21 1806 (1980). 
[13] J.M. Kosterlitz and D. Thouless, J. Phys. C 6, 1181 

(1973); J.M. Kosterlitz, J. Phys. C 7, 1046 (1974). 
[14] D. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 1201 

(1977). 

[15] In this formulation it is not necessary to follow the com- 
plex motion and reconnection of the pairs with each 
other, as such motion is reflected statistically in the 
change of the distribution function after each time step. 
Our recombination process at a is the same as the "bulk 
recombination" process discussed by other authors who 
follow the collision of individual pairs RLZ]. 

[16] H. Yano, T. Jocha, and N. Wada, Phys. Rev. B 60 543 
(1999). 

[17] N. Antunes, L. Bettencourt, and W. Zurek, Phys. Rev. 

Lett. 82, 2824 (1999). 
[18] T. Shaw et al, Phys. Rev. Lett. 76, 2551 (1995). 



5 



[19] C. Muzny and N. Clark, Phys. Rev. Lett. 68, 804 (1992); 
N. Clark, private communication. 



